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ABSTRACT 


Aims. A complete census of obscured Active Galactic Nuclei (AGN) is necessary to reveal the history of the super massive black hole 
(SMBH) growth and galaxy evolution in the Universe given the complex feedback processes and the fact that much of this growth 
occurs in an obscured phase. In this context, hard X-ray surveys and dedicated follow-up observations represent a unique tool for 
selecting highly absorbed AGN and for characterizing the obscuring matter surrounding the SMBH. Here we focus on the absorption 
and reflection occurring in highly luminous, quasar-like AGN, to study the relation between the geometry of the absorbing matter 
and the AGN nature (e.g. X-ray, optical, and radio properties), and to help to determine the column density dependency on the AGN 
luminosity. 

Methods. The Swift/BAT nine-month survey observed 153 AGN, ah with ultra-hard X-ray BAT fluxes in excess of KT 11 erg cm -2 
s -1 and an average redshift of 0.03. Among them, four of the most luminous BAT AGN (44.73 < LogL B AT < 45.31) were selected 
as targets of Suzaku follow-up observations: J2246.0+3941 (3C 452), J0407.4+0339 (3C 105), J0318.7+6828, and J0918.5+0425. 
The column density, scattered/reflected emission, the properties of the Fe K line, and a possible variability are fully analyzed. For 
the latter, the spectral properties from Chandra, XMM-Newton and Swift/XRT public observations were compared with the present 
Suzaku analysis, adding an original spectral analysis when non was available from the literature. 

Results. Of our sample, 3C 452 is the only certain Compton-thick AGN candidate because of i) the high absorption (Nn ~ 4 x 10 23 
cm -2 ) and strong Compton reflection; ii) the lack of variability; iii) the ’’buried” nature, i.e. the low scattering fraction (< 0.5%) and 
the extremely low relative [OIII] luminosity. In contrast 3C 105 is not reflection-dominated, despite the comparable column density, 
X-ray luminosity and radio morphology, but shows a strong long-term variability in flux and scattering fraction, consistent with the 
soft emission being scattered from a distant region (e.g., the narrow emission line region). The sample presents high (> 100) X-to- 
[OIII] luminosity ratios, with an extreme value of R^ in j ~ 800 for 3C 452, confirming the [OIII] luminosity to be affected by residual 
extinction in presence of mild absorption, especially for ’’buried” AGN such as 3C 452. Three of our targets are powerful FRII radio 
galaxies, which is shown by their high luminosity and absorption; this makes them the most luminous and absorbed AGN of the BAT 
Seyfert survey despite the inversely proportional N H - L x relation. 

Key words. X-ray: Galaxies - Galaxies: Active 


1. Introduction 

Absorption by circumnuclear gas and dust surrounding the ac- 
creting supermassive black hole (SMBH) is a common feature 
of active galactic nuclei (AGN). According to the unified model 
(Antonucci 1993), different observational features can be ex- 
plained on the basis of the AGN orientation with respect to the 
observer, such that an object is classified as unobscured Type I 
AGN if the central broad line region appears in the optical spec- 
trum, while the lack of these features occurs if the nuclear re- 
gions are seen at larger inclination angles and the line of sight 
intercepts optically thick cold matter distributed in a toroidal ge- 
ometry. This model is supported by both X-ray (Bassani et al. 
1999, Tozzi et al. 2006, Cappi et al. 2006) and infrared (IR) ob- 
servations (Alonso-Herrero et al. 2006, Fiore et al. 2008), where 
a high percentage of AGN have X-ray column densities greater 
than 3 x 10 22 cm -2 . 

If the bulk of the SMBH growth is expected during an obscured 
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phase of the AGN evolution (Hopkins et al. 2006), and its accre- 
tion is tightly connected to the galaxy bulge evolution and star- 
burst activity in a feedback process (Ferrarese & Merritt 2000, 
Brandt & Hasinger 2005), the complete census of the obscured 
AGN is a fundamental key to unlock the SMBH growth his- 
tory. In the local Universe, a significant fraction of Seyfert 2 
(Risaliti et al. 1999, Guainazzi et al. 2005) are Compton-thick 
(Nr > 1.5 x 10 24 cm -2 ). Current synthesis models of the cos- 
mic X-ray background (CXB), given by the contribution of un- 
obscured and obscured AGN integrated over the cosmic time 
(Setti & Woltjer 1989, Comastri et al. 1995), predict that a high 
percentage of Compton-thick AGN contribute to the energy den- 
sity peak at 30 keV (Gilli et al. 2001). 

Despite the high penetrating power of X-ray photons, the na- 
ture of this population is still elusive because of the strong 
obstacles against their detection, because for column densities 
above 3 x 10 23 cm -2 the X-ray flux from the nucleus is re- 
duced by more than a factor of 10. In this context, hard X- 
ray surveys above 10-15 keV performed in recent years by 
Swift/BAT (14-195 keV, Tuelleretal. 2010) and INTEGRAF 
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(10-100 keV, Beckmann et al. 2009) represent a unique tool for 
obtaining a complete, unbiased sample of AGN for column den- 
sities < 3 x 10 25 cm -2 . The need for X-ray selected AGN surveys 
is also of particular interest if we consider that more than 50% of 
the X-ray AGN are not detected at other wavelengths and prob- 
ably trace different AGN populations and/or evolutionary stages 
(Hickox et al. 2009). The discovery of a new class of “buried” 
AGN (Ueda et al. 2007) from the BAT hard X-ray survey is a 
direct consequence of the new observing window. 

The full characterization of the obscuring matter surrounding 
the SMBH requires more than the Nh evaluation. Other impor- 
tant observational clues are found along the entire X-ray spec- 
trum: from the large equivalent widths of the Fe K line to the 
flattened slope and the broad bump peaking around 20-30 keV 
due to Compton down- scattering of high-energy photons. These 
additional diagnostics for obscuration are necessary to avoid 
an underprediction of the column density that could affect the 
computed X-ray luminosity even in the hard E> 10 keV band 
(LaMassa et al. 2011, Melendez et al. 2008). 

The Hard X-ray Detector (HXD) on board the NASA/JAXA 
Suzaku mission extends the bandpass of the X-ray focusing tele- 
scope up to 600 keV with an unprecedent sensitivity in the 10- 
100 keV energy range thanks to the low in-orbit background 
(Takahashi et al. 2007). 

The Swift/BAT nine-month survey observed 153 AGN, all with 
ultra-hard X-ray BAT fluxes (14-195 keV) in excess of 10“ 11 
erg cm -2 s -1 and an average redshift of 0.03. In this paper we 
analyze the follow-up Suzaku observations of four of the most 
luminous BAT AGN (44.73 < Log L BAT < 45.31). 

The reason for the choice to focus on absorption and reflection of 
high X-ray luminosity AGN is the dependency of the obscured 
fraction on the AGN luminosity. The difficulty in mapping the 
space density of highly obscured AGN in the local Universe is 
overcome by the large uncertainties in determining the column 
density distribution as a function of luminosity and cosmic time, 
which is strictly correlated to the accretion history and evolu- 
tion of AGN. While there is still debate on the redshift distribu- 
tion of obscured AGN (La Franca et al. 2005, Ballantyne et al. 
2006, Treister et al. 2010), in the past decade X-ray (Ueda et al. 
2003, La Franca et al. 2005) optical (Simpson 2005) and in- 
frared (Fiore et al. 2009) studies seem to confirm the inverse pro- 
portionality of the obscured population as a function of the lumi- 
nosity, with a moderately obscured to unobscured sources ratio 
decreasing from 4 to 1 forL x < 10 42 andL x > 1 0 45 , respectively 
(Gilli et al. 2007). In this context, characterizing the obscuring 
matter surrounding high-luminosity AGN helps to separate the 
various interpretations (Hasinger 2008, Diamond- Stanic et al. 
2009) and to better constrain the Nh-L x correlation. 

In this paper, the extended Suzaku bandpass is exploited to study 
within the same observation the soft/hard X-ray emission of 
four AGN from the nine-month Swift/BAT survey (Tueller et al. 
2008). The column density, scattered/reflected emission, and the 
properties of the Fe K line are fully analyzed along with their 
variability, while the nine-month-averaged Swift/BAT spectrum 
helps to better constrain the spectral shape in the hard X-ray do- 
main. 

2. Observations 

The targets of our analysis are Swift J2246. 0+3941 (3C 452), 
J0407.4+0339 (3C 105), J0318.7+6828, and J0918.5+0425. We 
summarize their basic parameters in Table 1 and also list the 
host galaxy, AGN radio, and optical classification collected 
from the literature. The first two targets, J2246. 0+3941 and 


J0407.4+0339, are well-known 3C catalog FRII Radio Galaxies, 
and we refer to them as 3C452 and3C 105, while J03 18.7+6828 
is an FRII Giant Radio Galaxy (GRG) with a linear size of 2 
Mpc. Optical spectra, available for all the sources, identify our 
targets as Type II AGN. The last, and farthest, J09 18.5 +0425 
has been identified as an AGN for the first time by the nine- 
month Swift/BAT survey, but was then not detected in the 4. So- 
cut of the 22-month survey, most likely because of variability 
(Tueller et al. 2010). 

The Japanese Suzaku (ASTRO-E2) X-ray mission, launched 
in 2005, achieves its wide (0.2-600 keV) X-ray energy range 
through combining four imaging CCD cameras (the X-ray imag- 
ing spectrometers, XIS) with a non-imaging, collimated hard X- 
ray detector (HXD). Each XIS detector, three front- (XIS-0,2,3) 
and one back-illuminated (XIS-1), is located in the focal plane of 
a dedicated X-ray telescope with a bandpass in the range 0.2-12 
keV, while the HXD consists of Si PIN photo-diodes and GSO 
scintillation counters, covering the 10-70 keV and 40-600 keV 
energy range (Mitsuda et al. 2007). 

Suzaku observed our sources (in HXD nominal pointing) be- 
tween June 2007 and February 2008 (see observation details in 
Table 2). Since the FI (front-illuminated) XIS-2 detector became 
inoperable in 2006 November, no XIS-2 data are available. We 
present the analysis of XISs and HXD/PIN data because the flux 
of our targets was too faint to be detected by HXD/GSO. In ad- 
dition to the Suzaku data, we also integrate the 14-195 keV BAT 
spectra in the spectral analysis, which were time-averaged over 
the nine-month survey to better constrain the fit. 

Except for 3C 105, which is characterized by shorter detections, 
the net (dead-time corrected) exposure spans from about 60 to 
80 ks for the XIS-BI, and from about 50 to 65 ks for the high- 
energy detector. 

3. Data analysis 

Observation products (spectra, images, light curves) were ex- 
tracted from the cleaned version 2.3.12.25 processed event files 
distributed by the Suzaku team, with the most recent HEAsoft 
package. The XIS source spectra and light curves, generated by 
combining the 3 x 3 and 5x5 edit modes in xselect, were 
extracted from circular regions with a radius ranging from 2.2' 
to 4.2'. The minimum 2.2' radius region was used for 3C 452, 
allowing us to take more than 80% of the point source photons 
(see Suzaku technical description) while avoiding contamination 
from external sources in the field. A 4.2' radius region, which 
collects 99% of the point source flux, was used for the other 
sources because there are no apparent confusing sources (see 
Sec. 4.2) inside the extracting region. Background spectra were 
extracted from nearby free-emission regions of about 2.2'-5' ra- 
dius. After producing the rmf and arf response matrices with the 
xisrmfgen and xissimarfgen tools provided by the Suzaku 
team, we combined the two XIS-FI spectra and responses with 
addascaspec for faster spectral fits, grouping (with grppha) 
both XIS-FI and XIS-BI source spectra using a factor greater 
than 25 cts/bin. 

Because HXD/PIN is a non-imaging detector, a model of the 
time-variable non X-ray background (NXB) is provided by the 
HXD team (dead-time corrected) with a systematic uncertainty 
of about 1.3% at lcr confidence level in the 1 5^40 keV band 
(for a 10 ks exposure, see Mizuno et al. (2008)). Since the NXB 
event rate is ten times higher than the real background, the par- 
ticle background spectra and light curves were increased by a 
factor of 10 with the fparkey and fcalc tools to suppress 
the Poisson error. A common good time interval (GTI) between 
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Table 1. 


List of Targets 0 


Source 

R.A. 

Dec. 

Redshift 

Log L B at^ 

Optical type 6 

Radio type J 

Host Galaxy 6 

Ref/ 

3C452 

341.4532 

39.6877 

0.0811 

44.73 

Sy 2 

FRII RG 

E 

1,2,3 

3C 105 

61.8186 

3.7071 

0.089 

44.83 

Sy 2 

FRII RG 

? 

1,2,3 

J03 18.7+6828 

49.5791 

68.4921 

0.0901 

44.85 

Syl.9 

FRII GRG 

? 

4, 5 

J09 18.5+0425 

139.5011 

4.4184 

0.1564 

45.31 

QSO 2 

/ 

Merg? 

4 


° The position and redshift (J2000 coordinates) are taken from the NASA/IPAC Extragalactic Database. 
b Logarithm of the luminosity in the 14—195 keV band in erg s _1 from the BAT nine-month AGN survey. 
c FRII RG = Fanaroff-Riley Class II Radio Galaxy; GRG = Giant Radio Galaxy. 

d GRG = Giant Radio Galaxy; Sy = Seyfert; QSO = Quasar; / = No radio counterpart is associated to the source. 

6 The AGN host galaxy is indicated as E (elliptical), ? (optical image too faint to be classified), Merg? (possible galaxy merging). 
f 1 = Laing et al. (1994), 2 = Jackson & Rawlings (1997), 3 = Black et al. (1992), 4 = Tueller et al. (2008), 5 = Schoenmakers et al. (1998). 

Table 2. 


Suzaku Observation Log 


Source 

Obs ID 

Date 

Exposure 0 

Ct. Rate 0 

3C 452 

702073010 

2007-06-16 

66696, 57383 

0.037, 0.039 

3C 105 

702074010 

2008-02-05 

38319, 45518 

0.031,0.024 

J03 18.7+6828 

702075010 

2007-09-22 

77140, 65631 

0.091,0.027 

J09 18.5+0425 

702076010 

2007-11-04 

61019, 51899 

0.024, 0.006 


° Exposure time (s) and background-subtracted count rate (cts s 1 ) for XIS-BI (back-illuminated) and HXD/PIN respectively in the 0.5-10 keV 
and 10-30 keV band 


the PIN cleaned event file and the NXB was created to obtain 
dead-time corrected spectra and light curves. Like the NXB, the 
PIN also still suffers from CXB contamination (about 5%). This 
was estimated from the PIN response file for the flat emission 
distribution (provided by the HXD team) and then simulated 
with xspec using the HEAO-1 CXB spectrum (Boldt 1987). 
The particle background and CXB spectra were combined us- 
ing mathpha. The PIN response files, available from the Suzaku 
CALDB as generated in June 2009, were chosen according to 
the HXD nominal pointing and the observation epoch. 

4. Spectral analysis 

Using the X-ray spectral fitting package Xspec vl2 (Amaud 
1996), three main models were tested to describe the primary 
emission from the AGN. Each target was initially fitted with the 
simplest case; if the model was unable to describe the spectrum 
or the fit was significantly improved by introducing additional 
components, we adopted more complex models. In addition to 
the continuum emission, a Gaussian profile was added to char- 
acterize the Fe Kcr fluorescence line at 6.4 keV (zgauss) as the 
result of reprocessing in the surrounding material (e.g., accre- 
tion disk, torus, ionized gas). The physical width was kept free 
to vary. If the spectrum showed no significant fluorescence line, 
we derived an upper limit of the equivalent width, fixing the line 
centroid (in the source rest frame) at 6.4 keV. 

We used the zphabs code to model the photoelectric absorp- 


tion by cold matter across the line of sight with cross sec- 
tions from Balucinska-Church & McCammon (1992). The neu- 
tral absorption taking place in our galaxy was accounted for 
(even when not explicitly mentioned) with the tbabs model 
(Wilms et al. 2000). The galactic column density, from the HI 
map of Dickey & Lockman (1990), was provided by the nh pro- 
gram of the HEAsoft package and is listed in Table 3. Solar 
abundances as given by Anders & Grevesse (1989) were as- 
sumed throughout the analysis, with cosmological parameters of 
H 0 = 70 km s -1 Mpc -1 , £2 M = 0.27 and Q vac = 0.73. 


4.1. Simple model 

The X-ray continuum emission from a Type II AGN is character- 
ized, as a first-order approximation, by an absorbed power law 
(Mushotzky et al. 1993), resulting from thermal Comptonization 
of soft-seed accretion-disk photons by a hot diffuse corona 
(Haardt & Maraschi 1991). Defining as N H and <x(E) the column 
density and photoelectric cross section respectively of the cold 
matter surrounding the black hole, this simple model takes the 
form F(E) = exp{-Nncr(E)} E _r (T represents the photon index), 
translated into zphabs*zpowerlaw in Xspec terminology. 

We assumed the absorber to be neutral, although ionized ab- 
sorbers in radio loud AGN have recently been found (e.g., in 
the BLRG 3C 445, Braito et al. (2011)), because the data quality 
does not allow such a detailed analysis. 
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4.2. Partial-covering model 

Several Type II AGN spectra show, in addition to the absorbed 
primary emission, a secondary less absorbed, or not at all ab- 
sorbed, emission at lower energies (Turner etal. 1997). Many 
interpretations could be applied to this observational feature: 
it could be 1) the result of contamination by other sources 
(e.g., X-ray binaries/starburst emission, Maiolino et al. (1998)); 
or 2) a clumpy, dusty absorber that partially covers the nucleus 
(Malaguti et al. 1999); or 3) AGN light scattered back into our 
line of sight by hot gas (Matt etal. 1996); or 4) the sum of 
unresolved emission lines from photoionized gas (Bianchi et al. 
2006). 

Given the point spread function (PSF) of the XIS instruments, 
we cannot rule out other sources on the basis of the Suzaku 
imaging capabilities alone. For this reason, both the Swift/XRT 
and, when available, the XMM-Newton and Chandra images of 
our targets were analyzed to investigate contamination by point 
sources or the lobes of the associated radio galaxy. 

If we interpret the secondary component as intrinsic emission 
from the AGN, this complex spectrum can be described by the 
partial-covering model (see Holt etal. 1980), where a fraction 
of the primary emission escapes without being absorbed. It re- 
quires the hydrogen column density and the covering fraction 
of the partial absorber (f c ), and can be written as F(E) = {f c 
expj-NHcr(E)} + (1 - f c )} E -r . This model is useful to describe 
not only a clumpy absorber but also scattered nuclear emission 
escaping from different regions, characterized by lower column 
density, and redirected to our line of sight, assuming the same 
slope for the scattered power law. We refer to the value 1 - f c 
as the scattering fraction f scatt of the absorbed power law. The 
partial-covering code takes the form zpcfabs*zpowerlaw. If 
f c < 1, a fraction (1 - f c ) of the primary emission escapes with- 
out being absorbed, while a f c = 1 returns the simple absorbed 
model. 

The soft X-ray emission could also be thermal, due to collision- 
ally ionized diffuse gas that could be present in the host galaxy 
or nearby the AGN itself (Isobe et al. 2002), or the sum of unre- 
solved lines produced by photoionized gas (Reeves et al. 2010). 
The XIS spectral resolution does not allow us to separate the two 
sources, and we only searched for additional thermal gas with 
the apec component to better constrain the scattered/partially 
covered power law. 

4.3. Double power law model 

In general, the second, scattered power-law could be absorbed 
by less dense clouds and be characterized by a different slope 
than the direct emission. The third model increases the com- 
plexity of the fit by adopting two different absorbed power-laws 
of the form F(E) = exp{-N^cx(E)} E“ r ‘ + exp{-N^cr(E)) E“ r \ 
where a and b mark distinct values, or zphabs*zpowerlaw + 
zphabs*zpowerlaw in Xspec. 

4.4. Compton reflection 

Given the extended energy range of the PIN, the additional 
contribution of Compton reflection by optically thick matter 
could arise in the hard part of the spectrum. The Xspec pexrav 
code models the angle-dependent Compton reflection of an inci- 
dent exponentially cut-off power law by a plane-parallel, semi- 
infinite medium of cold electrons that approximate an accre- 
tion disk (Magdziarz & Zdziarski 1995). We applied this code 
to a more general case where the AGN primary radiation could 


be reflected by optically thick matter present in the accretion 
disk or in the surrounding torus/absorber (see the caveats in 
Murphy & Yaqoob (2009)), resulting as the sum of the incident 
cut-off power-law and the reflection component. The amount of 
reflection depends on both geometry and composition of the re- 
flector, parameterized in the model by the inclination of the re- 
flecting material and the gas abundance, with the possibility of 
varying the Fe abundance alone. The model also returns the rel- 
ative strength R (rel re fi in the code) of the reflection component 
to that of the incident cut-off power law component, defined as 
R = Cl/ 2/r, where Cl is the solid angle covered by the reflector. If 
R = 0, the pexrav model translates into a simple exponentially 
cut-off power law. Each model (simple model, partial-covering 
model, double power-law model) was tested with and without 
the zphabs*pexrav code for all sources. When a significant 
reflection component was not found from the spectral fitting, 
we evaluated an upper limit to the relative reflection by fitting 
the spectra with a partially covered pexrav model of the form 
zpcfabs*pexrav. As previously noted, for R = 0 this model 
turns into a partial-covering model. 

5. 3C 452 

The source 3C 452 is an FRII radio galaxy with a symmetri- 
cal double-lobe morphology and a total angular extent of about 
l'x4 / (Black et al. 1992). Because it does not lie in a rich cluster, 
it is not contaminated by thermal X-ray emission related to the 
intracluster medium (Isobe et al. 2002, hereafter 102). From the 
Hubble observation, no UV excess (Wills et al. 2002) and no op- 
tical central compact core are detected (Chiaberge et al. 2002), 
while the NICMOS near-IR instrument reveals a faint compact 
source 1” southwest of the nucleus (Madrid et al. 2006), con- 
firmed by the weak mid-IR core detected by Spitzer (Ogle et al. 
2006). 

The Chandra/ACIS instrument observed the source on August 
2001: a detailed spectral analysis of the two lobes and central 
AGN emission is reported by 102, and the results for the AGN 
properties are confirmed by Evans et al. (2006). From 102, the 
soft X-ray emission from the extended radio structure is fitted by 
a power law, interpreted as inverse compton (IC) and significant 
in the 2-5 keV band, with the addition of the Raymond-Smith 
(RS) thermal emission in the 0.5 - 2 keV band from hot, diffuse 
gas associated to the host galaxy. Hodges-Kluck et al. (2010) 
confirmed this emission throughout the region covered by the 
radio lobes, although with the non-thermal emission dominant 
near the hot spots, but they associated the thermal emission to 
the radio galaxy cocoon or a bounding shock and not to the ac- 
tual intergalactic medium. 

5. 1. Extracting the AGN emission 

In the analysis of the Suzaku/XIS image, the required extracting 
region of the central point source covers part of the extended ra- 
dio galaxy because the PSF is wider than that of Chandra/ACIS: 
the AGN emission is contaminated by the diffuse emission at 
energies lower than 5 keV. This can be clearly seen in the top 
panels of Figure 1, which shows the Chandra/ACIS and XIS1 
0.5-2 keV smoothed image of 3 C 452 in the left and right pan- 
els, while the white circle refers to the 2.2' radius extracting re- 
gion. The 1 .4 GHz Very Large Array (VLA) radio map (Laing 
R.A. unpublished) is superimposed with green contours: the soft, 
cocoon shaped X-ray emission traces the radio contours in the 
ACIS image while it is not visible in the 2-10 keV energy selec- 
tion. 
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Fig. 1. Chandra/ ACIS (left panel) and Suzaku/XISl (right panel) 
images of 3C 452 in the 0.5-2 keV (top panels) and 2-10 keV 
(bottom panels) energy ranges. The 1 .4 GHz VLA radio map is 
superimposed with green contours and the white circle indicates 
the Suzaku 22' radius extracting region. 


Because we aim to analyze the AGN spectrum, it is necessary to 
separate the X-ray photons of the central core from the extended 
emission. The lobe-induced contamination in the XIS spectrum 
was checked by simulating the diffuse and core detection of the 
XIS1 on the basis of the ACIS observation with the xissim 
tool. We extracted the ACIS total spectrum from the same re- 
gion on which the Suzaku 3C 452 analysis is based: the central 
core counts were taken from a 6” circular region, while the ex- 
tended emission was obtained from subtracting the core region 
from the total one (Figure 2, top panel). The two Chandra spec- 
tra were fitted according to the models and best-fit parameters 
reported in 102, leaving the normalization free to vary. Applying 
the XIS1 response matrix to the resulting AGN and extended 
emission spectral models, we obtained a count rate of a fake 
XIS 1 detection of the two components. 

The bottom panel of Figure 2 shows the simulated XIS1 de- 
tection (the AGN emission in red and the extended emis- 
sion in green) and, as a comparison, the real observation (in 
black), which is the sum of the diffuse and central core emis- 
sion. The simulation reproduces the real detection well, i.e., the 
Suzaku/XIS spectra are contaminated by the diffuse emission be- 
low 3 keV. To model the AGN emission, two methods are possi- 
ble: 1) the XIS spectrum is analyzed only above 3 keV, losing all 
information about the possible scattered/partially covered com- 
ponent; 2) the model of the soft X-ray emission, based on the 
102 results, is added to the AGN model, and the full energy band 
XIS spectrum (0.5-10 keV) is analyzed. We selected this last 
option to be able to also explore the possible partially-covered 
emission. 

5.2. Spectral fitting 

We analyzed the spectrum in the full energy band by adding the 
XIS (FI and BI in the 0.5-10 keV energy range), the PIN (10- 


Black = Total, Red = AGN, Green = Extended emission 




Energy (keV) 

Fig. 2. Top panel: Chandra/ ACIS spectrum of the 2.2' radius 
region centered on the AGN (black line), the AGN alone (red 
line), and the extended region after subtracting the AGN core 
(green line). Bottom panel: the Suzaku/XISl observed spectrum 
is shown in black, while the red and green lines refer to the XIS1 
simulated spectra of the AGN and the extended emission. 


30 keV), and the nine-month-averaged BAT (14-195 keV) data. 
The normalization between XIS and PIN flux was accounted 
for using a 1 . 1 8 factor in the PIN model, given by the cross- 
correlation of the Crab observation in HXD nominal pointing, 
while the BAT normalization was kept free. Because the ex- 
tended emission is not expected to vary during the six year time 
interval between the Chandra and Suzaku observations, we can 
model the extended emission on the basis of the Chandra spectral 
parameters: zpowerlaw + raymond (the thermal emission fit- 
ting is based on the calculations of Raymond & Smith (1977)), 
with a power- law photon index of 1.68 and a gas temperature 
of 1.36 keV, as reported in 102. The relative flux for the ex- 
tended emission derived from spectral fitting between ACIS and 
Suzaku agrees within the cross-calibration error of the two in- 
struments. The simple absorbed power law leaves an excess in 
the PIN spectrum, and the reflection component is needed. In 
the pexrav model we could not determine the high-energy cut- 
off, which was fixed at 200 keV, or the inclination and chemical 
composition of the reflector, fixed at 63° (the default value) and 
solar abundance. The partial-covering model did not improve the 
fit, given the model-crowded soft band, but it allowed to fix an 
upper limit to the scattering fraction, which is below 0.5%. A po- 
tential absorption of the secondary power law could not be con- 
strained and was removed from the model. The best-fit model is 
zpcfabs*zpowerlaw + zgauss + pexrav, and the main pa- 
rameters are reported in Table 3. 
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3C 452 



Energy (keV) 


nated by the southern hotspot and the wealth of other knots and 
structures suggests (Leahy et al. 1997) that the system is rapidly 
evolving, with a short lifetime. Chandra detected X-ray emis- 
sion from the core and southern hot spot of 3C 105 in December 
2007 (Massaro et al. 2010, hereafter M10), both from the region 
where the jet appears to enter the hot spot and at the terminal 
hotspot itself. No extended, thermal X-ray emission is detected. 
A faint, unresolved near-IR source is visible 4" to the west- 
southwest off the nucleus (Madrid et al. 2006), while we were 
unable to find any information on the UV emission from the lit- 
erature. Because the exposure is shorter than that of the other 
sources, we extracted the counts from a 4.2' radius region that 
also includes the southern hotspot, to collect up to 99% of the 
source photons. 



1 10 100 
Energy (keV) 


Fig. 3. Folded (top panel) and unfolded (bottom panel) spectra 
and best-fit models of 3C 452: black = XIS-FI, red = XIS-BI, 
green = HXD/PIN, blue = Swift/BAT. 


5.3. Fitting results 

The source 3C 452 is highly absorbed, with a local column den- 
sity N h of 4 x 10 23 cm" 2 . The hard energy band of the spectrum 
is dominated by Compton reflection of the intrinsic power-law 
photons onto a cold, thick matter, and the source seems to see 
the absorber covering a solid angle greater than 4/r: the relative 
strength of the reflection component to the transmitted one, R, is 
higher than 400. This value is unphysical, but it can be mathe- 
matically explained as the result of an inhomogeneous absorber, 
with a significant fraction of the solid angle covered by a gas 
thicker than that along the line of sight (Risaliti et al. 2002). The 
nuclear emission is flattened, with a power-law photon index of 
1.55, while the average Seyfert II slope ranges from 1.7 to 1.9 
(Turner & Pounds 1989, Nandra & Pounds 1994). The relative 
normalization of the BAT spectrum is close to 1 : the lack of vari- 
ability is consistent with 3C 452 behaving as a Compton- thick 
AGN (see Section 9). However, the Fe Kcr line, centered at 6.43 
keV (in the source rest frame) and narrow, shows an equivalent 
width (EW) of only 164 eV, in contrast with what is expected 
from Compton-thick sources (see Sec. 12). 


6. 3C 105 

The source 3C 105 is an FRII narrow-line radio galaxy with 
a spatial size of the radio emission of 764 kpc (Baum et al. 
1988, Hardcastle et al. 1998), associated to an elliptical host 
galaxy (Donzelli et al. 2007). From the radio map of Leahy et al. 
(1997), both jets are detected, but the radio emission is domi- 



3C 105 

Black = XIS-FI, Red = XIS-BI, Green = PIN, Blue = BAT 



1 10 100 
Energy (keV) 


Fig. 4. Folded (top panel) and unfolded (bottom panel) spectra 
and best-fit models of 3C 105: black = XIS-FI, red = XIS-BI, 
green = HXD/PIN, blue = Swift/BAT. 


6.1. Spectral fitting 

Following the procedure reported in Sec. 4, the X-ray spectrum 
of 3C 105, given by the sum of the XIS, PIN, and BAT data, ini- 
tially fitted by a simple absorbed power-law model, which leaves 
a strong residual in the soft energy band. The addition of ther- 
mal emission (with the apec model) requires a gas temperature 
that is too high (> 2 keV) to fit the data, while a typical value 
of 0.65 keV is observed in Seyfert 2 galaxies (Guainazzi et al. 
2005). The best fit is given by the partial-covering model, with 
the scattered power-law modeling the excess at softer energies. 
Using a double power-law, with different column density and 
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Table 3. 


Suzaku Spectral Fit Compared Results 


Source 

n H a 

[10 22 cm -2 ] 

N h 6 

[10 22 cm -2 ] 

r 

f c 

^scatt 

Fe K E d 
[keV] 

Fe K EW rf 
[eV] 

R e 

BAT 
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A 2 /d.o.f 
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0.12 

43.52!$f 
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1 + 
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< 0.5% 
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no 

164*45 
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p p 

116.76/175 

3C 105 

0.12 

45.96*3 

O 0\ 

o o 
+ 1 

OO 

L- 

1-4*05% 

6 40 +0 - 07 
u, ^-.0.05 

04 

ON 
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< 1.8 

o o 
+ 1 

OO 

NO 

68.18/97 

J0318. 7+6828 

0.35 

5-26t« 

OO OO 
p p 
© © 
+ 1 

IT) 

IT) 

i-2 :°oi% 

6.4(f) 

/tq +32 

63-27 

< 1.4 

o 

00 

1 + 
p p 

167.68/172 

J0918. 5+0425 

0.035 

16.33 2 _fi 3 

i n^+0.20 

/z -0.19 

0.9!";®% 

6.4(f) 

< 98 

<2.6 

i 74 +1 - 1 6 
A * /h -0.87 

79.21/132 


a Galactic column density provided by the nh tool, using the HI map by Dickey & Lockman (1990). 
b Source intrinsic column density. 

c Fraction of the scattered component relative to the intrinsic power law. 

d Centroid energy and equivalent width respect to the whole continuum of the Fe fluorescence K a line at the rest frame of the source redshift. 
e The relative strength of the reflection component to the transmitted one, defined as R= O/ In, where Q is the solid angle of the reflector viewed 
from the nucleus. 

f Normalization factor of the BAT flux respect to the XIS detection. 


photon index, the same set of spectral parameters was found. 
The best-fit model is zpcfabs*zpowerlaw + zgauss (Fig. 4). 
While the addition of Compton reflection does not improve the 
fit, we can evaluate an upper limit to the relative reflection. We 
fit the spectrum with a partially-covered pexrav model of the 
form zpcfabs*pexrav + zgauss. The composition and line 
of sight angle are settled following the model default values (so- 
lar abundance composition and 63° inclination angle). If R is 
equal to 0, the pexrav model translates into a partially covered 
power law, which is our best-fit model. Fixing the column den- 
sity, scattering fraction and normalization to the best fit values, 
we calculated the relative reflection confidence range as a func- 
tion of the power law slope (see Table 3, Col. 7). 

6.2. Fitting results 

Adding the nine-month-averaged BAT spectrum to the Suzaku 
observation, a normalization factor of 1.7 for the BAT spectrum 
is needed, i.e., the BAT flux is about 60% higher than the Suzaku 
detection in the same band, which implies a possible long-term 
variability (see Sec. 9). The source 3C 105 intrinsic power-law 
emission is highly absorbed by an Nr ~ 5 x 10 23 cm -2 , while 
1 .4% of the nuclear X-ray photons is scattered at lower energies 
and reaches the observer unabsorbed. 

As previously noted, the XIS extracted region also covers the ra- 
dio galaxy southern hotspot. According to M10, its X-ray emis- 
sion is much more fainter than the core: its flux is only about 
30% of the core in the 1-2 keV energy range and lower than 
2% in the 2-7 keV band. However, because we interpret the soft 
X-ray excess as scattered nuclear emission, we checked for a 
possible external contamination by comparing their fluxes. Table 
4 reports a summary for all sources of the flux and luminosity 
emitted in the main bands with (left value) and without (right 
value) absorption (galactic and local). In the 0.5-2 keV energy 
range we obtain an absorbed flux of 50 x 10 -15 erg cm -2 s _1 . 
From M10, the total flux of the southern hotspot region in the 
0.5-2 keV energy range is 5 x 10 -15 erg cm -2 s -1 , ten times 
lower than the Suzaku/XIS observation. If there is contamina- 
tion, it falls within the uncertainty of the scattering fraction. 

The spectral shape of 3C 105 differs from the XIS analysis of 


Winter et al. (2009a) (W09a hereafter), especially in the power 
law slope and scattering fraction. This is due to the addition 
of the PIN data, which allows us to constrain the photon in- 
dex to a steeper value. Although the source is not dominated 
by Compton reflection in the hard band, we obtained a reflection 
relative strength upper limit of 1.8, translating into a covering 
solid angle less than 3.6n. A narrow 136 eV EW Fe K a line 
centered at 6.4 keV is found. 


7. Swift J031 8.7+6828 

The source J03 18.7+6828 is identified as an FRII GRG (see 
Schoenmakers et al. (1998) for a full description of the radio 
and optical properties). The optical spectrum (INT telescope in 
La Palma) is typical of a narrow line AGN, while the optical 
image shows a faint galaxy-like object associated to the radio 
core but the resolution is not enough to clarify the host mor- 
phology. We could not find any information on the IR or UV 
maps from literature, even if the IR emission, mostly in the K s 
band, is found within 30" from the source position in the Two 
Micron All Sky Survey (2MASS) at IPAC. XMM-Newton ob- 
served J03 18.7+6828 on February 2007 (a detailed analysis can 
be found in Winter et al. (2008), hereafter W08): the source im- 
age, from the XMM/pn event file, in the 0.5-2 keV energy range 
does not show the presence of external sources and/or contami- 
nation by thermal emission. With a radio galaxy angular size of 
15' (2 Mpc at z = 0.09), the giant radio lobes are well outside 
the Suzaku/XIS extracting region. 

7.1. Spectral fitting 

A simple, absorbed power law can fit the source J03 18.7+6828 
with sufficient accuracy, even in the soft 0.5-2 keV energy band, 
but the best-fit was gained by modeling its emission with a par- 
tially covered power law. A Gaussian was then added to model 
the faint Fe fluorescence line, fixing its rest energy at 6.4 keV, 
and the best-fit model takes the form zpcfabs*zpowerlaw + 
zgauss (Fig. 5). No significant reflection is found in the hard 
part of the spectrum: we evaluated the upper limit to a possible 
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Table 4. 


Flux and Luminosity 0 




Flux in 10“ 

-12 -2 -1 

erg cm s 


Log Luminosity in erg s 1 

Source 

0.5-2 keV 

2-10 keV 

10^10 keV 

14-195 keV 

0.5-2 keV 

2-10 keV 

10^10 keV 

14-195 keV 

3C452 

0.01, 1.79 

1.85,6.59 

14.46, 14.97 

37.78, 38.01 

41.20, 43.45 

43.47, 43.88 

44.36, 44.38 

44.78, 44.78 

3C 105 

0.05,5.11 

2.11, 8.28 

9.00, 9.82 

22.92, 23.26 

41.99, 44.00 

43.62, 44.21 

44.25, 44.29 

44.66, 44.66 

J03 18.7+6828 

0.06, 2.55 

4.47, 6.03 

10.23, 10.34 

32.19, 32.24 

42.08, 43.71 

43.96, 44.09 

44.32, 44.32 

44.81,44.81 

J0918. 5+0425 

0.01, 1.45 

1.55,2.87 

4.05,4.15 

11.33, 11.37 

41.84, 44.00 

44.03, 44.29 

44.44, 44.45 

44.89, 44.89 


a Fluxes and luminosity are reported with both the galactic and local absorption (left value of each column) and as intrinsic, not absorbed, 
emission (right value of each column) for the following energy bands: 0.5-2 keV, 2-10 keV, 10-40 keV and 14-195 keV. 




1 10 100 
Energy (keV) 


Fig. 5. Folded (top panel) and unfolded (bottom panel) spectra 
and best-fit models of Swift J03 18.7+6828: black = XIS-FI, red 
= XIS-BI, green = HXD/PIN, blue = Swift/BAT. 


Compton hump following the same procedure applied to 3C 105 
(see Sec. 6.1). 

7.2. Fitting results 

The nuclear emission is absorbed by a column density of Nr ~ 
5 x 10 22 cm -2 , about an order of magnitude lower than the pre- 
vious 3C sources, and the partial covering model returns a scat- 
tering factor of 1.2%. The power-law photon index appears to 
be flatter than the mean value associated to Seyfert galaxies as 


3C452. If reflection is present, the spectral fitting reports an up- 
per value of 1.4 to its relative strength, i.e., the absorbing matter 
subtends a solid angle lower than 2.87T. 


8. Swift J091 8.5+0425 

The source Swift J09 18.5+0425, the farthest among our sample 
(z = 0.15644) and the only radio-quiet AGN, can be defined as 
peculiar for many reasons. Detected by Rosat/PSPC in March 
2000, the source has been identified as AGN only after the nine- 
month Swift/BAT survey. After being classified as a QSO Type 
II (Tueller et al. 2008) AGN, J0918. 5+0425 has been removed 
from the second BAT survey release because of variability, and 
it is not part of the 58-month catalog. The Suzaku and Swift X- 
ray spectra are the only X-ray data available, at present, while 
no additional studies in the other bands are published, except 
for the catalog of Swift/BAT AGN optical spectra (Winter et al. 
2010, see Sec. 10). Its optical counterpart from the SDSS 1 shows 
a peculiar double-lobed galaxy-like object, possibly the site of a 
galaxy merging, but the resolution is not high enough to separate 
its components. 


8.1. Spectral fitting 

The PIN extracted spectrum results are not consistent with either 
the XIS or BAT spectra. This could be related to a low statistical 
significance due to the low count rate or to an incorrect parti- 
cle background model, which contributes about 95% of the total 
count rate. To investigate the latter hypothesis, we computed the 
ratio between the total detected count rate and the NXB: because 
the source emission is strongly affected by the particle back- 
ground, we should obtain a constant value. 

As shown in Figure 7 (top panel), after 9 x 10 4 s (spacecraft 
time), the ratio deviates from the constant behavior. When we 
remove the corrupted time interval, the PIN exposure time de- 
creases from 52 to 43 ks, with an average count rate of 0.01 cts 
s" 1 . The ratio resulting from the cleaned PIN light curve (Figure 
7, bottom panel) can now be well fitted (x 2 /dof = 18.10/18) by 
a constant function, and the PIN spectrum is consistent with the 
XIS and BAT data (Figure 6). Following the same procedure as 
described for the other sources, the best-fit model is given by a 
partially covered power-law of the form zpcfabs*zpowerlaw. 
Adding a Gaussian profile to model the Fe K a emission, we can 


1 Sloan Digital Sky Survey (SDSS), www.sdss.org 
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1 10 100 
Energy (keV) 


Fig. 6. Folded (top panel) and unfolded (bottom panel) spectra 
and best-fit models of Swift J09 18.5 +0425: black = XIS-FI, red 
= XIS-BI, green = HXD/PIN, blue = Swift/BAT. 



Time (s) 


Fig. 7. Background-unsubtracted PIN (at the top), background 
(in the middle) light curves and their ratio ( at the bottom) of 
J09 18.5+0425 are plotted before (top panel) and after (bottom 
panel) removing of the background-corrupted time interval. The 
goodness of the ratio fit by a constant function is also reported 
for the two cases. 


evaluate an upper limit to the line intensity, fixing the centroid 
rest frame energy at 6.4 keV. 


8.2. Fitting results 

The X-ray spectrum of J09 18.5+0425 was modeled by a par- 
tially covered power law, resulting in a 0.9% scattered factor. 
The source is moderately absorbed (Nh ~ 2 x 10 23 cm -2 lo- 
cal column density) and shows a weak Fe fluorescence line (EW 
< 90 eV). We stress that this is the highest luminosity source in 
our sample with an intrinsic, unabsorbed 2-10 keV luminosity 
of 2 x 10 44 erg s -1 , confirming the QSO nature of the source as 
reported by the Swift/BAT AGN catalog. 

9. Variability 

Strong variability is found in the X-ray emission of AGN, from 
intrinsic changes in the source flux to variability in the spectral 
properties (column density, photon index) and absorbed/emitted 
features (Mushotzky et al. 1993). The study of these variances 
and their time scales can provide information on the size 
and geometry of the emitting X-ray region (McHardy et al. 
2006), the physical conditions of the X-ray reprocessing media 
(Miniutti et al. 2007) and/or the physical properties of the ab- 
sorbing material close to the black hole (Miller et al. 2008). 

In contrast, if the line of sight to the source is blocked by ob- 
scuring matter, the direct emission is absorbed and only the pho- 
tons scattered by distant, cold matter can reach the observer, los- 


ing any trace of variability. A Compton-thick AGN, where lit- 
tle to no high direct emission escapes below 10 keV while the 
hard energy part of the spectrum is dominated by reflection, can 
be traced by the lack of variability, both in flux and spectral 
shape. This is the case of Circinus, with no significant change 
during nine years of observations (Winter et al. 2009b and refer- 
ences therein), and the mildly Compton-thick source NGC 4945 
(Itoh et al. 2008). 

There are two types of variability that we can test on the ba- 
sis of the Suzaku observation: variability on the time scales 
of the observation (ks) and variability between different ob- 
servations (time scales of months/years). In the first case, we 
examined the light curves of both PIN and XIS detections to 
look for a significant variation in the source count rate, while 
long-term variability was studied by comparing the spectral 
parameters and fluxes reported here with other observations 
(see Sobolewska & Papadakis (2009) for a review). Except for 
J09 18.5+0425, we also checked for a possible variability in 
the BAT 14-195 keV count rate along the 66-month survey 
(Baumgartner et al. 2011): the discovery of variability in the en- 
ergy range not dominated by reflection, i.e., the hard band ob- 
served by the BAT, in contrast with the lack of flux changes in 
the tens of keV, would be an additional evidence for the presence 
of Comp ton- thick material. 
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9.1. Short-term variability 

The XIS light curves, plotted with a 1 ks and 5.76 ks (the orbital 
time scale) binning, were extracted in two energy ranges, 0.5-2 
keV (Figure 8) and 2-10 keV (Figure 9). This allowed to explore 
the variability of the soft X-ray emission. 

For 3C 452, XIS light curves were only extracted in the 3-10 
keV band because of contamination in the soft, extended emis- 
sion. The PIN is plotted with a 5.76 ks binning along with the 
NXB light curve because it is strongly affected by the flaring be- 
havior of the particle background. We do not show the PIN light 
curve for J09 18.5+0425 since the NXB contributes the 97.6% 
of the observation and a longer exposure is needed to evalu- 
ate a possible variability. All XIS light curves were fitted by a 
constant function using the QDP fitting suite. No variability was 
found in the soft band but, given the low level of counts, this re- 
sult cannot be unambiguously related to the reprocessed nature 
of the emission. Above 2 keV, we found a slight deviation from 
a constant trend for the sources 3C 105 and J03 18.7+6828 in the 
5.76 ks binning, where we also obtained the lowest upper limit 
for reflection. We confirm that no strong, significant short-term 
variability (factor of 2 or higher) is found below 1 0 keV. The 
PIN light curves instead show a remarkable count rate temporal 
variation, which, however, reproduces the NXB fluctuations: the 
temporal behavior of the PIN is clearly affected by the particle 
background, hiding any potential AGN variability. 

9.2. Long-term variability 

The spectral properties obtained from the analysis of Chandra, 
XMM-Newton, and Swift/XRT data were collected from the lit- 
erature (see Table 5). For 3C 105, a simple model to which 
we applied, differences in the photon index or column density 
could fall within the uncertainty of our best-fit model. The X- 
ray spectral analysis for other sources has not been published 
yet (e.g., the XMM-Newton observation of 3C 452). While a 
detailed analysis of these observations is not the aim of this 
work, we compared our results with those from Suzaku by fit- 
ting these spectra with our best-fit model, letting the normal- 
ization free to vary for changes in luminosity and adding the 
nine-month-averaged BAT data to constrain the photon index. 
When the source could not be described by our best-fit parame- 
ters, we looked for intrinsic variations in the spectral parameters, 
with the assumption that our best-fit model reflects the physi- 
cal properties of the AGN. In the case of the XRT data analy- 
sis the Fe K line was not modeled because of the low-energy 
resolution of the instrument. The Swift/XRT data filtering was 
performed by means of the HE AS ARC xrtpipeline follow- 
ing the standard event screening criteria using the CALDB re- 
lease of the observation epoch. The Chandra/ ACIS analysis is 
based on the cleaned event files of the Chandra archive third re- 
processing (ASCDSVER = 7.6.11.3, CALDB = 3.4.2). For the 
XMM-Newton/pn data, we reduced the data following the stan- 
dard procedures as adviced by the XMM-Newton SOC using 
SAS v. 10.0.0. The extracting regions range from 5b" -60" for 
the XRT images to 30" and 1" for XMM-Newton and Chandra 
images, respectively. 

9.2.1. 3C 452 

The spectral parameters of 3C 452 agree well with the de- 
tailed analysis of the Chandra observation performed by 102 
and Evans et al. (2006). They also confirm the reflection com- 
ponent and the narrow Fe K feature. We included the thermal 


emission as modeled by 102 (accounted for by the normaliza- 
tion of the extended emission in the Suzaku analysis) and the 
scattering component in the analysis of the XMM-Newton ob- 
servation (the Swift/XRT detection was not available). A par- 
tially covered power law coupled with strong reflection gives a 
good fit (^ 2 /d.o.f = 296/291) to the data for the XIS analysis. 
We found a low fraction of scattered emission (< 0.9%) with an 
uncertainty due to the presence of thermal photons. The spectral 
properties of 3C 452, i.e., high column density, flattened photon 
index, narrow Fe K line, and strong reflection, do not show vari- 
ability during six years. This results, coupled with an unchanged 
absorption-corrected 2-10 keV luminosity, confirm 3C 452 as a 
highly absorbed, reflection-dominated source. 

9.2.2. 3C 105 

The X-ray analysis of 3C 105 has been published by 
Ajello et al. (2008) (Swift/XRT + BAT) and Massaro et al. 
(2010) (Swift/XRT, Chandra/ACIS and XMM-Newton/pn). 
While these authors generally confirm the strong absorption 
(> 3 x 10 23 cm -2 ), the use of simple absorbed power laws to 
model the spectrum and the lack of information regarding the 
source luminosity did not allow them to correctly evaluate a 
possible intrinsic variability. After summing the four available 
Swift/XRT observations, which span a time interval of five days, 
we try to describe the Swift/XRT, Chandra/ACIS and XMM- 
Newton spectra with the best-fit model and parameters obtained 
from the Suzaku spectral analysis (see Sec. 6.2). The source is 
intrinsically variable, and the new best-fit parameters are listed in 
Table 5 for each observation. The lack of reflection and the nar- 
row Fe K line are confirmed but the scattering fraction is lower 
than 0.2% for Chandra and Swift observations, and lower than 
0.5% for the XMM-Newton detection. These changes are cou- 
pled to an absorption-corrected 2-10 keV luminosity that de- 
creases by more than 50% from the middle of 2006 to the first 
months of 2008. A slight variability (~ 10% - 20%) in the col- 
umn density is also found. The intrinsic variable nature of 3C 


3C105 

Black = Swift/XRT, Red = Chandra/ACIS, Green = Suzaku/XIS, Blue = XMM/pn 



Fig. 11. Unfolded spectrum of 3C 105 as observed by Swift/XRT 
(2006-07-11/16, black line), Chandra/ACIS (2007-12-17, red 
line), Suzaku/XIS (2008-02-05, green line), and XMM-Newton 
(2008-02-25, blue line). 


105 is clearly visible in Fig. 11, where the unfolded spectra of 
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Table 5. 


Spectral Fit Results from Chandra/XMM-Newton/Swift XRT 
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\ 12 m 

1 ‘ ' -0.19 

3.3^- 9 % 

6.4(f), 0.044!q 04 

N 

1.5(obs) 

3 

J0318. 7+6828 

SX 

2006-03-29 

A 

/ 

/ 

N 

0 . 8 (obs) 

3 


SX 

S+B 

2006- 04-05 

2007- 09-22 

A 

C 

52 6 6 M 
5 - 26 -0.41 

/ 

i.2!°!;% 

/ 

6.4(f), 0.063 

N 

N 

l.O(obs) 

1 . 2 (unabs) 

3 

4 


a Observatory name: S = Suzaku/XIS, B = Swift/BAT, C = Chandra/ACIS, X = XMM-Newton/EPIC, SX = Swift/XRT. 
b The best-fit models: A = N H *PL, B = N H *PL + Gauss, C = Pcfabs*PL + Gauss, D = Pcfabs*PL + Gauss + Pexrav, E = N H *PL + Gauss + 
Pexrav + Th, F = Pcfabs*PL + Gauss + Pexrav + Th. All models include the galactic absorption, not shown here for simplicity. 
c Source-intrinsic column density and primary power-law photon index (p indicates a value pegged at the lower/upper limit of 1.5/1. 9). 
d Where modeled, the fraction of the scattered component relative to the intrinsic power law. 

e Centroid energy and equivalent width with respect to the whole continuum of the Fe fluorescence K a line at the rest frame of the source 
redshift. 

f Y = present, N = not found, / = not modeled 

g (obs) = observed, (unabs) = absorption corrected, (pow) = power law. 
h 1 = Isobe et al. (2002), 2 = Evans et al. (2006), 3 = Winter et al. (2009a), 4 = here. 


the four observations are plotted together: a higher intrinsic lu- 
minosity is associated to lower values of the scattering fraction. 

9.2.3. Swift J031 8.7+6828 

A detailed analysis and comparison of the Swift/XRT and 
XMM-Newton observations of J03 18.7+6828 is reported in 
W08: flux variability is found, with the XMM detection showing 
the source to be brighter by ~ 30% two months before the two 
XRT observations. Suzaku observed the source more than one 
and a half years later: column density, power law slope and the 
Fe K a feature are consistent with the previous results, but the 
scattered emission at lower energies is significantly lower than 
the XMM-Newton observation. The observed luminosity in the 
2-10 keV resulting from our analysis is 0.9x 10 44 erg s _1 , consis- 
tent with the values given by the XRT observation. We confirm 
the long term variability reported by W08: the source is seen to 
be more luminous and more scattered in early 2006, followed 
by a decrease in X-ray power of about 30% after two months(no 
information is available for the scattering component in the XRT 
data), and the same luminosity is found more than one year later. 
The scattering factor shows a strong decrease from the XMM to 
the Suzaku observation. 


9.2.4. Swift J091 8.5+0425 

As previously noted, the Suzaku observation of J09 18.5 +0425 
represents the first attempt to study its broad band X-ray emis- 
sion after the discovery of its high hard X-ray luminosity in the 
BAT survey. Except for a detection in the soft band by Rosat, 


the Swift/XRT view J09 18.5+0425 is the only available observa- 
tion one can use to look for variability. Unfortunately, the count 
level is low (46 counts in the whole band) and it is not statisti- 
cally significant below 3 keV. We were able to fit the data with 
our best-fit model after removing the Gaussian profile. However, 
this is probably due to the uncertainty of the XRT detection and 
not to the lack of variability. Indeed, the flux in the 2-10 keV 
band is about 50% lower than the XIS observation and almost 
four times lower than the nine-month averaged BAT flux extrap- 
olated in the same energy range. We conclude that the source 
shows long-term variability, but longer exposures are needed for 
a detailed analysis. 

9.3. BAT 66-month lightcurves 

The Swift/BAT 66-month Crab weighted light curves in the 14- 
195 keV energy range of 3C 452, 3C 105, and J0318.7+6828 
(J091 8.5+0425 was removed from the BAT AGN catalog) were 
binned with the rebingausslc tool using a three-month bin, 
and are plotted in Figure 12. The horizontal continuous and 
dashed lines refer to the average rate and the +lcr standard devi- 
ation. Although the count rate of our sample, especially 3C 452, 
seems to increase/decrease in a coherent way, all points and the 
associated error bars lie within 2cr. 

10. Scattering fraction 

All sources show a secondary power law arising in the soft X-ray 
band. A prominent soft X-ray emission is a common feature of 
optically selected samples of Seyfert 2 galaxies, with a scattering 
fraction in the range 3%-10% (Turner et al. 1997, Cappi et al. 
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Table 6. 


Optical Properties of the Sample" 


Source 

Obs. Date 

E(B-V) g / 

EfB-V),,/ 

L(Hu) 

L(H/i) 

L[OIII]T5007 

L[OIII]/L(H/3) 

L[2-10 keV] C 

rX d 

^[OIII] 

3C452 

2006-11-21 

0.14 

0.84 

40.94 

40.09 

40.89 

6.31 

43.81 

832 

3C 105 

2006-11-21 

0.48 

0.92 

41.15 

40.26 

41.5 

17.38 

44.21 

513 

J0318. 7+6828 

2006-11-21 

0.72 

0 

41.17 

40.74 

41.64 

7.94 

44.08 

275 

J0918. 5+0425 

2004-03-09 

0.04 

0.26 

41.63 

41.03 

42.1 

14.75 

44.30 

159 


a De-reddened emission line properties taken from Winter et al. (2010), based on observations from the 2.1 m Kitt Peak National Observatory 
(KPNO) telescope except for J09 18.5+0425, which is based on archival SDSS data. All luminosities are expressed as logarithms in units of erg 
s" 1 . 

b Milky Way and source-intrinsic extinction. 

c Logarithm of the intrinsic 2-10 keV power law luminosity in erg s -1 . 

d The ratio of the intrinsic 2-10 keV power law luminosity with respect to the corrected L[OIII]. 


2006, Guainazzi et al. 2005). Follow-up Suzaku observations of 
Swift/BAT AGN led to the discovery of a new type of AGN 
(“buried” or “hidden” AGN) with an extremely low scattering 
fraction (0.5%, see Ueda et al. (2007), Eguchi et al. (2009)). A 
quarter of 32 AGN analyzed by Noguchi et al. (2010) from the 
XMM-Newton Serendipitous Source catalog presents scattering 
fractions below 0.5%. W09a found that hidden AGN constitute 
a high percentage of the BAT AGN sample (~ 24%), adopting 
the criteria that a hidden source is one where the scattering frac- 
tion (fscatt) is < 3% and the ratio of soft (0.5-2 keV) to hard 
(2-10 keV) observed flux is < 4%. This subsample of hidden 
AGN also includes 3C 105 and J091 8.5+0425, while 3C 452 
and J03 18.7+6828 are not part of this class because they have a 
fscatt >3%. 

We confirm the W09a findings for the first two sources, but on 
the basis of the Suzaku and XMM-Newton analysis, we infer an 
upper limit of 0.5% to the scattering fraction for 3C 452, prob- 
ably because W09 did not model the reflection component. The 
scattering fraction of J03 18.7+6828 is variable around the clas- 
sification limit. Considering that the W09 hidden AGN show an 
averaged column density Nr > 3 x 10 23 cm -2 , about an order 
of magnitude higher than the value shown by J03 1 8.7+6828, the 
object with the lowest column density in our sample, we do not 
classify J03 18.7+6828 as a “hidden” AGN. 

It is thought that the soft X-ray emission is most often due to 
scattered light by the photo-ionized gas in the opening part of 
the torus while the direct view to the nucleus is absorbed, and 
the relative intensity of the scattered component with respect to 
the intrinsic one can be related to the opening angle of the op- 
tically thick torus. This new class of objects (Ueda et al. 2007) 
is deeply buried with a small opening angle torus (< 20° for 
the most extreme cases) and viewed in a face-on geometry. As 
a consequence, the same object seen edge-on would appear as a 
heavily Compton-thick AGN (N H ~ 10 25 cm -2 ). Only 3C 452 
of our sample, which is strongly absorbed, reflection dominated 
and has a < 0.5% scattering fraction, could effectively fall within 
this new class of AGN. However, the three radio galaxies of our 
sample could be contaminated in the soft X-ray emission by the 
base of the jet in the unresolved nuclear region, so that our find- 
ings would only pose an upper limit to the scattering fraction 
(Hardcastle et al. 2009). 

Since the photo-ionized gas responsible for the scattered light is 
probably the same region that generates the [OIII] emission opti- 


cal lines (Bianchi et al. 2006) of the Narrow Line Region (NLR), 
a low torus opening angle could affect the optical band. This 
bias would explain the lack of hidden AGN in optically selected 
samples. Winter etal. (2010) (W10 hereafter) listed the optical 
properties of the BAT AGN sample, comparing the absorption- 
corrected [OIII]/15007 with the BAT 14-195 keV luminosity. 
Contrary to the results of Heckman et al. (2005), they found only 
a weak linear correlation between the reddening-corrected [OIII] 
luminosity and L x . They suggested that the [OIII] luminosity 
is affected by extinction poorly corrected by using the Balmer 
decrement, and thus it cannot be used as an unbiased and ac- 
curate tracer of the AGN power. Similar results were obtained 
by Noguchi etal. (2010) who found Rf 0llI \ > 10 (the ratio of 
the intrinsic 2-10 keV power-law luminosity with respect to the 
corrected L[OIII]) in their sample of < 3% scattering fraction 
AGN and a value >100 for the most extreme cases, which is 
higher than the 1-100 range that they found from a subsample 
of Seyfert2 from Bassani et al. (1999) selected in the same range 
of N h (0.6 - 20 x 10 23 cm -2 ). A strong test for the buried nature 
of our targets would be a high value of R^^ . 

In Table 6 the extinction-corrected luminosity of the main opti- 
cal emission lines of our sample is listed as given by W10, along 
with the Galactic and the intrinsic reddening applied, the latter as 
inferred from the Balmer decrement. The last column shows the 
ratio between the intrinsic 2-10 keV power law luminosity and 
the absorption-corrected [OIII] luminosity (Rj^ni])- All ratios are 
well above the 1-100 range of the Bassani et al. (1999) sample, 
and agree very well with the values of L[om]/L x for narrow line 
objects reported in W10. The scatter in R^^ in the sample is 
only a factor of 5, with 3C 452, the AGN with the lowest scat- 
tering fraction and a coherent hard X-ray luminosity, reaching 
the maximum with a value of R^^ ~ 800, which is close to the 
extreme cases of Noguchi et al. (2010). 

The clearest result from the sample optical properties is that the 
X-ray to [OIII] luminosity ratio follows the X-ray properties of 
our sources. The three radio galaxies present an increasing Rj^^ 
as the column density increases and the secondary, scattered 
emission at soft X-rays becomes dimmer. Although the sample 
is far from being complete, the [OIII] luminosity is still affected 
by extinction, according to the hypothesis of W10, and it is a 
biased tracer of the AGN power, with a bias level depending on 
the geometry of the surrounding matter. The lowest luminosity 
ratio indeed comes from the most peculiar object: the QSO-like, 
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Fig. 13. Top panel : The distribution in column density (N H ) and 
BAT 14-195 keV luminosity of the BAT nine-month survey 
Seyfert galaxies as given by Tueller et al. (2008). The red crosses 
and blue diamonds represent our sample column density from 
the BAT survey and the Suzaku analysis. Bottom panel. The 
BAT nine-month Seyferts subsample is defined by Log N H > 23 
and/or Log T\^—\^keV > 44.5 (dashed lines), i.e., the most lumi- 
nous and absorbed sources are selected. The symbols refer to the 
radio type as classified from the literature, the thicker symbols 
indicate the present sample. 


Fig. 12. Swift/BAT 66-month Crab weighted light curves of 3C 
452 (top panel), 3C 105 (central panel), and J03 18.7+6828 (bot- 
tom panel), plotted with a three-month binning. The horizontal 
continuous and dashed lines refer to the average rate and the ± 1 <x 
standard deviation. 


highly variable J091 8.5+0425, the only source that does not emit 
in the radio band and shows both the highest X-ray and [OIII] lu- 
minosity. 

Our findings seem to confirm the origin of the low 3C 452 [OIII] 
relative luminosity as induced by a clumpy, small opening an- 
gle absorber, blocking the emission of the NLR where the [OIII] 
lines arise, which proves the hidden nature of 3C 452. 


11. Radio connection 

As mentioned in the introduction, the amount of obscura- 
tion in AGN is inversely proportional to the X-ray luminos- 
ity (Gilli et al. 2007). Despite this observational evidence, our 
sources are highly luminous and highly absorbed. Here we ex- 
amine the possibility that their high column density could be 
related to their radio properties. 

In the sample, the most absorbed sources (3C 452 and 3C 105) 
are powerful FRII radio galaxies, J03 18.7+6828 is also a giant 
radio galaxy. The only source that is an exception from our as- 
sumption is J091 8.5+0425, which is a highly variable QSO re- 
moved from the BAT 9-month survey, however. 

We took into account the BAT nine-month catalog, and, after 
selecting the Seyfert galaxies, i.e., discarding the Blazar/BL- 
LAC/Galaxy sources, the column density as a function of the 
BAT 14-195 keV luminosity is plotted in Fig. 13 (top panel), 
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using the values reported by Tueller et al. (2008). The crosses 
refer to the BAT Seyfert catalog, the red crosses represent the 
present sample, while the blue diamonds show the column den- 
sity obtained from the Suzaku analysis. We assumed, as reported 
by W09a, a proportional relation between the BAT 14-195 keV 
and the softer Suzaku 2-10 keV luminosity. 

Evidently, our sources are at the edge of the Nr - L x relation, 

i.e., they are both the most absorbed and the most luminous of 
the BAT Seyfert catalog. If this peculiar feature is related to the 
radio power, we should not find sources similar in terms of radio 
properties in other regions of the plot. To investigate this possi- 
bility, we traced two arbitrary regions (see dashed lines in Fig. 
13) that define the most absorbed (Log Nr > 23) and/or the most 
luminous (Log Li 4 _i 9 5 k e v > 44.5) sources of the Seyfert sample. 
At the intersection we find our sources, except for J03 18.7+6828 
(the lowest diamond) which is, in fact, a Compton-thin AGN 
with a column density of Nr ~ 5 x 10 22 cm -2 . We searched the 
literature for the radio classification of the BAT sources selected 
in terms of high luminosity/column density, which was difficult 
because many sources are defined differently according to dif- 
ferent authors. The BAT Seyferts are divided into: i) radio quiet 
(including sources with a low radio flux that does not permit a 
radio classification); ii) radio jet (sources with a jet, but without 
a radio type/morphology definition); iii) broad line radio galaxy 
(BLRG); iv) flat spectrum radio galaxy (FSRG); v) giant radio 
galaxy (GRG); vi) Seyfert 2 FRII radio galaxy. The result is plot- 
ted in Fig. 13 (bottom panel) and the thicker symbols refer to our 
sample. 

It is not surprising to find out that different classes of radio galax- 
ies are grouped differently according to the X-ray luminosity and 
the amount of absorbing matter (see e.g. Sambruna et al. 1999), 
because the radio morphology depends on the inclination of the 
AGN with respect to the observer and is connected to the radio 
power. For example, the column density increases from BLRG 
(blue circles) to GRG (green triangles) and, finally, FRII Radio 
Galaxies (pink crosses). However, we note that, in addition to 3C 
452 and 3C 105, and considering the GRG as a different class of 
objects, there are only two other similar FRII galaxies in the se- 
lection, and both lie in the high-luminosity/absorption region. 
The only three GRG, including J03 18.7+6828, share the lumi- 
nosity range with the classic FRII, but they are all Compton-thin 
AGN. Radio-quiet sources, such as QSO J09 18.5+0425, are in- 
stead uniformly distributed in the Nr - L x plot. 

Given these findings, it is clear that our sources are four 
of the most luminous and obscured AGN of the BAT sur- 
vey because three of them are powerful FRII radio galaxies, 
while J091 8.5+0425 was confirmed in its peculiar nature. The 
Suzaku/XIS angular resolution does not allow us to separate the 
potential X-ray jet from the AGN core itself. However, if an X- 
ray emission enhancement were caused by the jet, the effect on 
the scattering fraction would be an overprediction of the parame- 
ter, which in turn would consolidate the hidden nature of 3C 452. 
A higher X-ray luminosity could also cause an increase of the X- 
ray/[OIII] luminosity ratio Rj^^ . While such an effect cannot be 
excluded, it must be noted that the highest Rj^^ is associated to 
the lowest 2-10 keV intrinsic luminosity (3C 452). 


12. Summary and discussion 

We studied four of the most luminous BAT AGN (44.73 < 
Log L BAT < 45.31): J2246.0+3941 (3C 452), J0407.4+0339 (3C 
105), J0318. 7+6828, and J0918.5+0425. Of these, 3C 452 is 


the only certain Compton-thick AGN candidate because of four 
properties: 

1. The source spectrum shows both high absorption (Nr ~ 
4 x 10 23 cm -2 ) and strong Compton reflection, which also 
flattens the intrinsic power-law photon index. The relative 
strength of the reflection is consistent with a complex ab- 
sorber where a large portion of the subtended solid angle is 
composed by gas that is optically thicker than the amount 
crossed along the line of sight; 

2. The XIS and PIN light curves show no significant short-term 
variability, and the spectral parameters of the source stay the 
same during six years. The lack of variability is a trace of 
reflection-dominated emission, because the nuclear intrin- 
sic variability is diluted by radiation coming from the torus, 
which covers a much more extended region. 

3. Although the soft X-ray band is contaminated by the ex- 
tended thermal gas and the lobes’ IC emission, a low scatter- 
ing fraction (< 0.5%) was found, confirmed by the XMM ob- 
servation for which the lower PSF size allowed us to extract 
only the AGN region. This aspect, coupled with an extremely 
high ratio of the intrinsic 2-10 keV to the [OIII] luminosity, 
classifies 3C 452 as hidden AGN, where, as reported in Sec. 
10, the source is seen face-on and the equatorial region is 
covered by Compton-thick matter. This scenario would be 
consistent with the hypothesis of a complex absorber. 

4. The EW of the fluorescence iron line is lower than ex- 
pected for a Compton-thick AGN and consistent with the 
evaluated column density, for an inclination angle > 60° 
(Ghisellini et al. 1994, Ikedaetal. 2009). However, such a 
low EW could be justified by a Fe abundance lower than 
the solar value, in particular, an EW = 164 eV is obtained 
for a Fe abundance in the range 0.2-0. 3 (Ballantyne et al. 
2002). The quality of the data only allowed us to evaluate 
an upper limit to the Fe abundance of 0.7. It should be noted 
that supersolar metal abundances are expected for broad line 
Seyfert Is (Komossa & Mathur 2001). 

These properties assume more importance when comparing 3C 
452 with 3C 105. Despite their comparable column density, red- 
shift, X-ray luminosity, and radio properties, 3C 105 shows a 
long-term variability and is not reflection-dominated, while 3C 
452 is. These differences underline the necessity of hard X-ray 
observations to better characterize the absorbing medium and to 
find sources. 

In contrast, the overall properties of 3C 105 lead us to classify 
this source as a Compton-thin AGN. The most prominent fea- 
ture is the strong variability in flux and scattering fraction found 
by comparing our results with the Chandra, XMM-Newton, and 
XRT observations. The flux variability is also confirmed by the 
high normalization factor between the PIN and BAT data (~ 1.7) 
and the slight variations in the Suzaku light curves. A luminosity 
variability, especially on time scales of months/years, is not sur- 
prising for Type II AGN not dominated by scattering, because 
the moderate column density in the line of sight only blocks the 
soft X-ray energy range (E < 3 keV). A change in the scat- 
tered/partially covered emission is instead very peculiar, and, 
coupled with the flux changes, it could reveal the distance and 
extension of the scattering region. Among the four observations, 
the highest scattering component (~ 1 .4%) was found in the XIS 
detection, which also needs a more extended extracting region. 
This could imply that we are including in the source model some 
external component. However, no external point sources or ex- 
tended thermal emission are detected from the analysis of the 
high-resolution Chandra/ ACIS image, while a possible contami- 
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nation of the southern hotspot would be below the scattering fac- 
tor value, considering also the statistical uncertainties (see Sec. 
6.2). In addition, the analysis of the XIS spectrum of 3C 105 re- 
ported in WO 8 returns a higher scattering fraction (~ 2.7%). If 
the scattering component is intrinsically variable, it can be inter- 
preted differently according to the proposed models: 

1 . The arising power law in the soft band is nuclear emission 
escaping through a partially covered medium. A change in 
its intensity can be associated to the motion of the clouds 
along the line of sight, according to the way Risaliti et al. 
(2002) modeled the variability in the absorbing column den- 
sity. Extending the computation of these authors to this case 
and using the 3C 105 black hole mass as given by W10, a 
ANr = 5x 10 23 and a time delay of 50 days (Chandra/Suzaku 
observations time lapse) leads to an absorber distance < 1 pc 
for a density p < 10 9 cm -3 . 

2. The soft emission is scattering from a distant region (e.g., the 
NLR). A change in the nuclear flux results in a change of the 
scattered emission, with a time delay given by the distance of 
the scattering region from the nucleus. An increase in lumi- 
nosity can also increment the ionization of the surrounding 
medium, increasing the scattering efficiency. 

Because we did not find absorption in the scattering component, 
it is unlikely that a fraction of the highly obscured intrinsic emis- 
sion can arise from the inner regions of the torus without absorp- 
tion. On the other hand, the second hypothesis correlates well 
with the high variability of the absorption- corrected X-ray lumi- 
nosity and the low L2-10 kev/[OIII] ratio. 

The same variability behavior is found in J03 18.7+6828, al- 
though weaker. This AGN is the least absorbed of our sample 
(Nh ~ 4 x 10 22 cm -2 ), and the only one not classified as a hid- 
den AGN. The scattering fraction is variable by around 3%, a 
change that correlates with a 30% variability in the 2-10 keV 
band. Compared with 3C 105, the scattering fraction from the 
XIS analysis is about half of the value from the XMM-Newton 
observation, although the first requires a wider extracting region. 
We conclude that the variability of the soft scattered emission 
seems a ubiquitous property of these two sources, but a dedi- 
cated X-ray monitoring for a longer time is needed to effectively 
relate these findings to the geometry and composition of the scat- 
tering/partially covering region. 

Finally, the farthest source of our sample, J09 18.5+0425, is a 
moderately absorbed Type II AGN. The strong variability that 
led this source to be removed from the 22-month BAT survey 
is confirmed from comparing the XRT, Suzaku and nine-month- 
averaged BAT data. However, longer exposures and a finer spa- 
tial resolution are needed for a more detailed analysis. 

Last but not least, some considerations on the absolute luminos- 
ity of the sample. The definition of a Type II QSO is still arbi- 
trary: according to Zakamska et al. (2003), the AGN requires a 
Log L[OIII]T5007 > 40.06 erg s -1 to be classified as QSO, while 
the QSO sample of Krumpe et al. (2008) is based on absorption- 
corrected, intrinsic 0.5-10 keV luminosities > 10 44 erg s -1 . The 
only source to satisfy both definitions is J09 18.5+0425, which 
is confirmed to be a Type II QSO as reported in Tueller et al. 
(2008). However, as pointed out in Sec. 10, the [OIII] luminos- 
ity could suffer from extinction in case of high absorbing column 
densities. If we only consider the X-ray classification, then 3C 
105 and J03 18.7+6828 can also be defined as QSO sources. 

The common properties of all sources, except J09 18.5+0425, are 
that they show QSO-like luminosity, are obscured, present low 
scattering fractions and are powerful FRII radio galaxies, the lat- 
ter playing a key role in the luminosity- absorption relation. Is 


this a secondary effect due to the unresolved jet base contam- 
inating the X-ray luminosity? Or is there an intrinsic, mutual 
connection between the j et production mechanism and the matter 
surrounding the SMBH? Given the limited number of sources, it 
is impossible to reveal the nature of this connection here, or if 
the decreasing number of Compton-thick objects at higher lu- 
minosities is related to the radio population at these energies. In 
this context, the increase in the AGN survey accuracy in terms of 
completeness and energy coverage, as proven by the Swift/BAT 
survey and follow-up optical and soft X-ray catalog, along with 
the improvement in the sensitivity and spectral resolution of hard 
X-ray observations (e.g. NuStar), will be a fundamental tool to 
unfold the AGN cosmic history. 
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Fig. 8. Suzaku/XIS light curves in the 0.5-2 keV energy range obtained with a 1 ks (top panels) and 5.76 ks (the orbital time, bottom 
panels) binning. The source 3C 452 is not shown given the thermal/lobes contamination in the soft band. 



Fig. 9. Suzaku/XIS light curves in the 2-10 keV energy range, except for 3C 452 cutted below 3 keV because of contamination, 
obtained with a 1 ks (top panels) and 5.76 ks (the orbital time, bottom panels) binning. 
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Fig. 10. Background-subtracted PIN andNXB light curves above 10 keV, except for J091 8.5+0425 where the 97.6% of the detection 
is affected by the NXB and a longer exposure is needed to evaluate a possible intensity variability. 


